Coulomb parameters and photoemission for the molecular metal TTF-TCNQ 
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We employ density-functional theory to calculate realistic parameters for an extended Hubbard 
model of the molecular metal TTF-TCNQ. Considering both intra- and intermolecular screening in 
the crystal, we find significant longer-range Coulomb interactions along the molecular stacks, as well 
as inter-stack coupling. We show that the long-range Coulomb term of the extended Hubbard model 
leads to a broadening of the spectral density, likely resolving the problems with the interpretation 
of photoemission experiments using a simple Hubbard model only. 
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As the first realization of a molecular metal, the quasi 
one-dimensional molecular crystal TTF-TCNQ has been 
studied thoroughly for more than thirty years P|. The 
different electro-negativity of the two molecules leads to 
charge transfer of about 0.6e from TTF to TCNQ, effec- 
tively doping molecular stacks of like molecules, which 
become conducting. Low dimensionality and strong 
Coulomb repulsion 0] lead to many-body effects which 
appear in spin susceptibility and angular-resolved 
photoemission experiments (ARPES) 4, 5|. Most promi- 
nently, TTF-TCNQ is one of the few systems in which 
spectroscopic signatures of spin-charge separation have 
been clearly observed . The interpretation of the mea- 
sured spectra relies, however, on rough estimates of the 
Coulomb parameters, which come from early experimen- 
tal 01 and theoretical 0, ^] works. They suggest that 
the Coulomb repulsion U between two electrons on a 
molecule is larger than the band-width W. Yet, de- 
spite their importance for a correct description of the 
electronic properties, more quantitative values for the 
Coulomb parameters of TTF-TCNQ are to date not 
available. 

In addition, a description of TTF-TCNQ in terms of 
a pure one-dimensional Hubbard model with only on- 
site Coulomb interaction seems to have problems. For 
the interpretation of recent ARPES data |^ the hop- 
ping matrix element t within a molecular chain has to 
be taken about twice of what is estimated from band 
structure 9] and experiment Moreover, using these 
parameters, the temperature dependence of the spectral 
function cannot be understood [lol . These calcula- 
tions, however, neglect the inter-molecular Coulomb re- 
pulsion y, the relevance of which was already suggested 
by Hubbard and through diffusive X-ray scattering 
experiments, which display Wigner-type fluctuations 

Our aim in this Letter is to provide a realistic model 
for the description of TTF-TCNQ and to understand the 




FIG. 1: (Color online) TTF-HOMO (left) and TCNQ-LUMO 
(right) of two neighboring molecules in TTF-TCNQ. Since 
both are planar 7r-orbitals, the molecules facing each other in 
a stack are quite close. This gives rise to the hopping t and 
also implies a strong nearest-neighbor Coulomb repulsion V. 

spectral function. To calculate the screened Coulomb pa- 
rameters for an extended Hubbard model, we proceed in 
two steps. Using density-functional theory (DFT), we 
first determine the screening due to the electrons on the 
same molecule (intra-molecular screening). To also in- 
clude the inter-molecular screening we describe the other 
molecules in the lattice by point polarizabilities and check 
this electrostatic approach against constrained-DFT cal- 
culations for small clusters. Our calculations yield a siz- 
able coupling between neighboring molecules on different 
stacks, as well as important longer-range Coulomb inter- 
actions V along the stacks. Focusing on the effect of V on 
the spectral density, we find that it leads to a broadening 
consistent with photoemission experiments. 

We start our calculation of the Coulomb matrix ele- 
ments from the bare Coulomb integrals between molec- 
ular orbitals (MOs) on two molecules a distance I apart 

VLre = Jd'rJd'r'poiv)piir')/\v-v'\, with V^,, = C/barc 

We calculate the charge density p{r) of the highest oc- 
cupied molecular orbital (HOMO) of TTF and the low- 
est unoccupied MO (LUMO) of TCNQ from all-electron 
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TTF TCNQ TTF TCNQ TTF TCNQ TTF TCNQ 

(7baro 5.9 52 3.1 2^9 K,,, 1.8 O y^l,e 1-2 1^2 

(7o 4.7 4.3 \/o 2.9 2.8 1.8 1.7 Vq" 1-2 1.2 

17 2 1.7 VI 0.9 V 0.55 0.4 V" 0.4 0.3 

TABLE I: Hubbard parameters for TTF-TCNQ. Uharc is the direct Coulomb integral, Uo includes the renormalization due to 
intra-molecular screening, and U is the screened on-site Coulomb term in the crystal. Correspondingly, V gives the first, V' 
the second, and V" the third nearest-neighbor Coulomb interaction between molecules on the same TTF or TCNQ stack. All 
energies are in eV. For comparison, the band-width is about 0.7 eV. 



density-functional theory (DFT) using Gaussian orbitals 
[l^ and the Perdew-Burke-Ernzcrhof functional [isj . 
The results for the on-site, J/baic, and the first, Vbaro, sec- 
ond, V^^j.^, and third, V^'^^.^, nearest neighbor Coulomb 
terms on the TTF/TCNQ stacks are listed in Table U 
We notice that Vbarc is substantial, being of the order 
of C/baic/2. This can be understood from Fig. ^ The 
TT-molecular orbitals stick out of the molecular planes so 
that the charges from the two molecules repel strongly. 
For distances I larger than twice the molecular radius, 
^'are givcn by c^/Z. As a side-product we obtain from 
the bonding-antibonding splitting for a pair of neigh- 
boring molecules the hopping matrix elements along the 
stacks: |i| = 0.15 eV for TTF and 0.18 eV for TCNQ. 

The on-site Coulomb parameter Uq including intra- 
molecular screening is obtained by calculating the DFT 
total energy of an isolated molecule, putting charge q 
in the HOMO or LUMO, respectively, and fitting to 
Eu{<l) =10 + + U^q^ /2. The results are listed in 
Table HI an d agree well with previous calculations for 
TCNQ [lj|. We note that intra-molecular screening re- 
duces J/baic by more than 1 eV. In order to obtain the 
Coulomb parameters Vq , we calculate the total energy of 
two molecules at distance Z, putting charge q/2 on each 
and fitting to Ev{q) = 2Eu{q/2) + bo + hq + Voq'^/A. In 
contrast to the on-site term, we find hardly any screen- 
ing of the Vbaro's, since the molecules mainly screen their 
own charge. 

Determination of the inter-molecular screening contri- 
bution requires calculations of the energy of an infinite 
lattice of molecules. Constraining a charge g to a single 
molecule allows to obtain U, while constraining it to two 
molecules at a distance I provides the different Vi . Unfor- 
tunately, such constrained DFT calculations can only 
be done for very small clusters of molecules. In order to 
reach the infinite-system limit, we therefore evaluate the 
reliability of approaches which represent the molecules 
by their polarizabilities. To determine the polarizability 
tensor of the isolated molecules, we calculate, with DFT, 
the dipole moments in a homogeneous external field along 
the principal axes and extract the linear response. The 
results are given in Table HTl and are consistent with esti- 
mates derived from the experimental dielectric constant 
through the Clausius-Mossotti relation. 

One might then replace a molecule in the lattice by 



f^xx ^yy O^zz 

TTF 226 160 88" 
TCNQ 440 184 82 

TABLE II: Principal elements of the polarizability tensor of 
the TTF and the TCNQ molecule in atomic units (bohr^). x 
(y) is the long (short) ajcis of the molecule, while z is perpen- 
dicular to the molecular plane. 

a point polarizability at the molecule's center of gravity. 
Such an approach, which works nicely, e.g., in the case of 
Cgo ll^, fails in the present case, since the molecules in 
a TTF/TCNQ stack come so close that they are inside 
the convergence radius of the dipole approximation. To 
avoid this problem we distribute the polarizability uni- 
formly over the non- hydrogen atoms of the molecule ■ 
To test the reliability of this approach, we compare the 
screening due to the nearest molecules, where the dis- 
tributed dipole approach should work least well, to con- 
strained DFT calculations, using an implementation of 
the method in DMol"^ 0|. Calculating the reduction SU 
of the on-site Coulomb term due to screening of just the 
two nearest neighbors of a TTF molecule in a stack, we 
find SU = 1 eV, while the distributed-dipole approach 
yields 0.9 eV. A similarly good agreement is found when 
taking also the second nearest neighbors on the stack 
into account: The constrained DFT calculation for this 5 
molecule system gives SU = 1.7 eV, compared to 1.6 eV 
from the distributed-dipole approach. This points at a 
quite small error, in particular when considering that the 
screening due to the more distant molecules along the 
stack will be increasingly well described by the dipole 
approximation. 

We therefore employ the distributed-dipole approach 
to compute the screening of the Coulomb interaction in a 
lattice of polarizable points with the experimental crystal 
structure of Ref. (IJj . We use the notation of Allen |2fl| , 
where a configuration of dipole moments on the lattice is 
denoted by a vector |p) and the dipole-dipole interaction 
on the lattice is described by a matrix T. If the dipoles 
arise via polarization, their energy in an external field 
l^cxt^ is given by {p\a-^ - T\p)/2 - {p\W^) , where a 
is the polarizability tensor. The configuration that min- 
imizes the energy follows from the variational principle: 
|p) = {a~^ — r)~^|E™*). As the energy is quadratic 
in the external field we can focus on pairs of electrons. 
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Writing the field of two point charges located at lattice 
sites n and m as E"''* = E„ + E^, we find that screening 
lowers the energy by 

SW = -(E„|(a-i - r)-i|E„) - (E„|(a-i - r)-i|E„} , 

where in the first term, which gives the screening of 
each individual charge, we have used translational in- 
variance, and in the second term, which describes the 
screening of the interaction, we have used inversion sym- 
metry. If the two point-charges are located at the same 
position, we get the screening of the on-site interaction 
SU = -(Eo|(a-i - r)-i|Eo), while for \n - m\ = I we 
find 6V^ = -(Eo|(a~i - r)-i|E,}. We have performed 
these calculations for clusters of up to 400 molecules. Ex- 
trapolating to the infinite lattice limit gives the results 
compiled in Table HI We find that for TTF-TCNQ the 
screening is very efficient, reducing U — Uq — SU hy more 
than 2 eV, while for the longer range interactions the 
screening is somewhat smaller. 

Comparing to the band-width W of about 0.7 eV 0, 
we find that U/W is larger than 2, and also the longer 
range interaction is substantial, with the nearest neigh- 
bor interaction V being about U/2. Even the interac- 
tion between third nearest neighbors is still compara- 
ble to W. The parameters for the TTF and TCNQ 
stacks are quite similar, with TTF, having smaller hop- 
ping and larger Coulomb parameters, being slightly more 
correlated. While the electronic structure in the ab- 
sence of inter-stack Coulomb interaction is essentially 
one-dimensional, we find that there is also a sizable cou- 
pling between neighboring molecules on different stacks. 
For two neighboring TTF stacks we find Vbarc = 1.6 eV, 
the same for neighboring TCNQ stacks, while the inter- 
action between TTF and TCNQ stacks is about 2.1 eV. 
Screening reduces all these values to about 0.4 eV, which 
is still comparable to the band-width. 

Our results therefore show that a realistic description 
of TTF-TCNQ requires an extended Hubbard model, in- 
cluding both inter-stack, as well as longer-range Coulomb 
interactions along the stacks. This is in stark contrast 
to existing pure one-dimensional Hubbard models with 
on-site term only {t-U-mode\) ^ J^, JA^. Due to the 
inter-chain coupling we expect Coulomb-drag effects , 
giving rise to interesting correlated dynamics of the indi- 
vidual stacks. Treating this coupling at mean-field level, 
one arrives at a one-dimensional Hubbard model, with 
Coulomb-parameters renormalized by the metallic neigh- 
boring chains. This will somewhat reduce U and V, while 
the longer-range interactions are suppressed more effi- 
ciently jjf- In the following we study the importance of 
the longer-range Coulomb interaction along the stacks by 
considering the extended Hubbard model of a single TTF 
or TCNQ chain 

H = -t^^ {c\^Cja + H.c.) + nii^Uii -f T^y^ n^nj 





^ ^^^^^^^ 
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FIG. 2: Spectral function for the one-dimensional t-U-V 
model on 20 sites. The full lines show the effect of V for 
a TCNQ chain with n = 0.6, t = 0.18, U = 1.7, and V in- 
creasing from to 0.9 eV. The dotted line is for a TTF chain 
with n = 1.4, t = -0.15, U = 2.0, and V = 1.0 eV. 



using the values for U and V given in Table {t-U-V- 
model). With filling n « 0.6 and 1.4 and very similar 
parameters, the models for a chain of TCNQ and TTF 
are closely related by an electron-hole transformation, 
with TTF being somewhat more correlated than TCNQ. 
The following discussion of the spectral function for a 
chain of TCNQ and applies therefore correspondingly to 
TTF. 

Since U is more than twice the band-width VF, it seems 
safe to assume that in a TCNQ-chain double occupancies 
are completely suppressed. The effect of V is then to 
maximize the distance between occupied sites. For filling 
0.6 this results in a Hubbard- Wigner state with occupa- 
tions . . . 1101011010. . . with a periodicity of 5 lattice sites 
7]. However, such a state would be an insulator with a 
gap opening with V . This is, indeed, what we find solving 
the t-U-V Hubbard chain for U ^ oo. For J7 = 1.7 eV, 
however, the chain is metallic, indicating that double oc- 
cupancies can not be neglected. Indeed, tor V — we 
find that the probability d of a doubly occupied site is 
d « 0.01, to be compared to d = 0.09 for U = 0. The 
effect of V is to further increase d. For V = 0.9 eV we 
have d « 0.027. To obtain the same value in the t-U 
model, we would have to double t from 0.18 eV to about 
0.37 eV. Evidently, the nearest-neighbor Coulomb term 
encourages hopping. To understand this, consider what 
happens when two electrons pass each other: first they 
must become neighbors (energy V), then, to pass, they 
must form a doubly occupied site (energy U), to finally 
end up in a configuration where the neighbors have ex- 
changed roles. This process thus involves effectively an 
energy U — V. This can be translated to an effective hop- 
ping tcs ~ tU I {U — V) for a corresponding tan-U model. 
For our parameters we find toS ~ 0.38 eV. 

Figure 121 shows the spectral function Aii{uj) of the t- 
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U-V chain as function of V. We see that the main effect 
of V is to broaden the spectrum around the Fermi level. 
A similar broadening has also been observed in Ref. 
and is reminiscent of what happens when we increase t 
in the t-U model as proposed in Ref. 0. Remarkably, 
the broadening can already be understood in first-order 
Rayleigh-Schrodinger perturbation theory: If we consider 
the V-tevia in the t-U-V Hamiltonian as a perturbation, 
then, to first order, the energy of the many-body states 
changes by ^ while the wave-functions re- 
main unchanged. From the Lehmann representation we 
thus see that the position of the spectral peaks shift lin- 
early in V. This shift is proportional to '^(^ij^in-iUj), 
which, for the low energy states, tends to increase with 
the energy. Thus the further a pole is from the Fermi 
level, the more it moves with V, leading to the observed 
broadening of the spectrum. 

We finally illustrate an interesting effect of the polar- 
izability a on the screening of the Coulomb interaction. 
For this, we consider a one-dimensional chain of polariz- 
able points with lattice spacing b, which we can solve an- 
alytically. For a < 6'^/4C(3), where C("-) is the Riemann 
zeta function, such a system is stable against spontaneous 
polarization, so we can use translational invariance of the 
dipole-dipole interaction matrix to diagonalize 1/a — T in 
a plane- wave basis (2flj . From the spectral representation 
we then obtain 



dV' = - 



27r6 



dk cos{lk) 



3(Li2(e*'=))2 



63 /4a - 5R(Li3(e*'=)) 



where Li2(z) and Li3(z) is the di- and trilogarithm, re- 
spectively |22|. From this solution we find that the 
screening is less efficient for larger I, even showing an- 
tiscreening for I > 2 and a not too close to the ferro- 
electric instability. Approaching the critical a, we find 
that the interactions become almost independent of 
the distance I. Thus close to the ferroelectric instabil- 
ity electronic correlations are lost. TTF-TCNQ should 
be very close to the transition. Hence it would be excit- 
ing to see what happens to the electronic correlations in 
TTF-TCNQ under sufficient hydrostatic pressure. 

In conclusion, we have calculated a realistic set of 
Coulomb parameters for TTF-TCNQ. We find that the 
commonly used values for the on-site interaction U are 
correct, but that interactions between electrons on neigh- 
boring sites are significant and must be taken into ac- 
count. Including the nearest neighbor interaction y on a 
TCNQ or TTF chain, we find that the spectra are broad- 
ened, an effect that can also be mimicked by increasing 
the hopping term t. Not having to assume an enhanced t 
to obtain the proper width of the low-temperature spec- 
tral function, should also lower the estimate of the tem- 
perature scale Tj, above which the signs of spin-charge 
separation are lost, to the experimentally observed range. 
In addition we have found that the Coulomb terms cou- 
ple different chains, possibly giving rise to interesting 



Coulomb-drag effects. 
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